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Abstract 

We formulate quantum optics to include frequency dependence in the modeling of optical networks. 
Entangled light pulses available for quantum cryptography are entangled not only in polarization but 
also, whether one wants it or not, in frequency We model effects of the frequency spectrum of faint 
polarization-entangled light pulses on detection statistics. For instance, we show how polarization en- 
tanglement combines with frequency entanglement in the variation of detection statistics with pulse 
energy. 

Attention is paid not only to single-photon light states but also to multi-photon states. These are 
needed (1) to analyze the dependence of statistics on energy and (2) to help in calibrating fiber couplers, 
lasers and other devices, even when their desired use is for the generation of single-photon light. 

PACS numbers: 03.65.-w, 03.65.Nk, 03.65.Ta, 84.30.Sk 
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PARTI 



1. INTRODUCTION 

The complexities of quantum optics, with its multiple integrals over frequency and wave 
vectors, tempt one to simplify, and indeed the groundbreaking equations that launched quan- 
tum key distribution (QKD) were simplified rather drastically, often leaving out altogether the 
frequency spectrum of the light involved. While on one hand the QKD equations have involved 
simplifications, on the other hand they invoke concepts of quantum decision theory, little used 
in quantum optics, such as trace distances between density operators as a measure of their dis- 
tinguishability. The motivation for putting the complications of frequency spectra back into the 
equations by which we model the faint light used in QKD comes from recognizing that both in 
implementing QKD systems and in designing eavesdropping attacks against them, frequency 
spectra play a crucial role. 

This report adapts quantum optics to deal as directly as possible with pulses of weak light 
propagating though optical fibers. The equations introduced here to model faint light give 
expression to frequency spectra, including frequency-entanglement; they also define and show 
examples of relevant partial traces of density operators for entangled, frequency-dependent 
light, needed to make use of an exceedingly useful relation between entangled- state QKD and 
QKD implemented without entanglement. 

The report grew from notes on techniques, some borrowed, others developed from scratch, 
needed to model a version of BB84 that uses polarization-entangled light. Polarization- 
entangled light from available sources is also frequency-entangled, and the driving question 
was how this frequency entanglement modulates the dependence of polarization-entangled 
QKD detection probabilities on mean photon number. 

Some subsequent papers dealing with frequency effects in polarization-entangled QKD, 
such as Ref. [1], use the techniques and results of this report, and, in particular use certain 
convolution integrals that are described in Sec. 10 and investigated in detail in Appendix C, 
with accompanying MATLAB programs given in Appendix F. 
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A. Quantum modeling 



By definition, quantum modeling invokes equations constrained in form to those of quantum 
mechanics, expressing a joint probability distribution of (theoretical) outcomes in terms of an 
initial density operator p at time t — 0, a hamiltonian evolution operator H, and a resolution 
of the identity consisting of a set of non-negative operators Mj satisfying J2j Mj < 1. These 
engender a probability of a (theoretical) outcome j [2] : 

Pr(j) = Tr[e- itH/h pe itH/h M j ]. (1.1) 

Whether or not one makes it explicit, at p, H, and Mj are functions of parameters that one 
views as under experimental control; all that a model can say is said in terms of how the prob- 
abilities depend, via p, H, and Mj, on these parameters. A system of equations for modeling 
particular devices includes equations that specify properties of the H, p, and Mj, thus special- 
izing the probability distribution Pr(j). The outcome j can be a list of components, e.g. one 
component for each of several detectors, in which case Pr(j) is viewed as a joint probability 
for the components of the outcome. 

In choosing quantum equations to model experiments with light, one expresses light by one 
or another density operator p, and one expresses detecting devices by operators Mj, possibly 
augmented by probe particles, as discussed in [2]. As is well known, the boundary between 
preparation of light and its detection is chosen by the modeler, and can be pushed around [3]. 
Different choices of p, H, or Mj set up different quantum models. Implicitly or explicitly, p 
is a function of variables that express the setting of various knobs on the laser and/or other 
devices that generate the light, and Mj is a function of variables that express knob settings 
on the detectors, such as those that control polarizing filters. To claim that a set of equations 
of quantum mechanics describes an arrangement of devices is to claim that the probabilities 
calculated from the equations more or less fit relative frequencies of experimental outcomes 
obtained, for some range of knob settings, from the arrangement of devices. 

As discussed in [2], choosing equations to model an experiment takes guesswork, indeed, 
two layers of guesswork. From experiments with devices and a first layer of guesswork [4, 5], 
one abstracts experimental relative frequencies that can be compared with modeled probabil- 
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ities. Demanding an approximate fit to experimental relative frequencies puts a constraint on 
probabilities of outcomes as functions of knob settings, and hence establishes a property of 
states and operators that can be judged as fitting an experimental situation. Still, diverse sys- 
tems of equations can always be found to agree with any given set of experimental relative 
frequencies and yet disagree among themselves in probabilities that they generate for arrange- 
ments of devices not yet explored. For this reason, arriving at equations of quantum mechanics 
by which to describe the measured behavior of QKD devices requires reaching beyond logic 
to make what may as well be called a guess. For this reason, the sensible use of equations in 
modeling QKD is hardly to 'prove security' but instead to help in achieving transmission of 
keys in the face of practical obstacles and to design eavesdropping attacks. 

By recognizing an irreducible freedom in choosing systems of equations to describe an ar- 
rangement of devices, we can clarify the mathematics of modeling, unencumbered by physical 
interpretations that are always subject to choice. Besides helping with QKD, the techniques 
of modeling presented here can also serve other uses of faint light that generate interesting 
joint detection statistics expressible in quantum mechanics but unknown to classical physics, 
correlations that allow the invention of new kinds of "cameras" with which to see and respond 
to the physical world. 

B. Aims in developing a framework 

Picture an experimental network involving light sources and detectors linked by fibers as 
shown in Fig. 1, where the blob in the middle can include phase shifters, fiber couplers, and 
attenuators, as well as conversions from fiber to free space and back again to fiber. Mainly I 
discuss so-called single-mode fibers; (most of these actually support two polarizations). The 
various fibers of a network need not be alike; for instance they can vary in their propaga- 
tion constants and in their attenuation. Although most of the discussion is in terms of fiber, 
free-space links can also be included. The mathematical framework offered allows modeling 
the variety of responses encountered experimentally in light detection [2]; in particular, the 
response to single-photon states need not be binary. 
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FIG. 1 : Optical network. 



The aim in developing this framework for the analysis of optical networks is this: 

1 . Provide equations to express violations of Bell inequalities. 

2. Provide for modeling networks assembled from smaller pieces, like tinker toys, by splic- 
ing fibers of one to fibers of another. 

3. Provide for convenient expression of pieces of networks that are free of frequency con- 
version but that have diverse fibers with diverse propagation constants at any single fre- 
quency. I.e provide mathematics convenient for expressing single-frequency modes for 
networks (or parts of networks) containing fibers that differ in their propagation con- 
stants. 



C. Approach 

We split up the task of modeling a quantum network into the following modules: 

1. Develop mathematics for the quantum mechanics of a set of uncoupled, lossless trans- 
mission lines, each line expressed by a set of modes, where each mode supports a range 
of frequencies propagating in two directions, denoted "+" and "— ". 

(a) Corresponding to each mode, introduce a creation operator, and define single-mode 
quantum states as the vacuum state acted on by superpositions of these operators 
over some frequency band. 
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FIG. 2: Dissipation modeled by coupling of desired modes to loss modes. 

(b) Extend to many modes. 

(c) Define projection operators, each of which corresponds to a subspace of n photon 
states on a Hilbert space that is a tensor product of some number of modes. 

i. Use lattice of projections to define subspaces. 

(d) Define detection operators on such states. 

(e) Perform calculations by use of commutation relations. 

2. Use scattering theory to approximate interacting transmission lines by operators that 
convert an in-state on non-interacting lines to an out-state on non-interacting lines. (See 



3. Model loss and dispersion by introducing coupling of desired modes to extraneous 
modes, as illustrated in Fig. 2. 

4. Model frequency preserving interactions among transmission lines by unitary transfor- 
mations acting on creation operators for all desired and extraneous modes. 

5. Sketch two applications to a network for quantum key distribution, in which these con- 
cepts and techniques work together. 

2. MODES, COMMUTATION RULES, AND LIGHT STATES 

In analogy with classical electromagnetics, we will model a path, such as an optical fiber, as 
a system of modes, with each mode supporting a continuous range of frequencies. We assume 



Fig. 1.) 
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a vacuum state |0), normalized so that 

(0|0) = 1. (2.1) 

Other states are defined by various creation operators acting on the vacuum state. These in 
turn are defined in terms of more singular single-frequency creation operators, as follows. For 
modes a, b, ... , one could introduce annihilation operators a, b, . . . ; however we avoid the 
clutter of the "hats," so that, for example, b expresses a mode name in some contexts but in 
other contexts expresses an annihilation operator for that mode. Each mode is bi-directional; 
we call one direction "+" and the other "— ". Let a\_(uj) be the creation operator for excitation 
at an angular frequency uj in the "+" direction, and let a_(oj) be the creation operator for 
propagation in the "— " direction. The operators a + and a_ are defined for non-negative uj. It 
is often convenient to combine the "+" and "— " operators into a single operator. Since these 
are both defined for uo > 0, we can define a(uo) for — oo < uo < oo by 

Ja+M forcu>0, 
a(uo) = < (2.2) 

I a_(— oo) for uj < 0. 
For any given mode, a single-photon state is any state of the form 

<4|0>, (2.3) 

where the creation operator aj. is defined by 

/oo 
duo f(uo)a) {uj) , (2.4) 
-oo 

and / is any square-integrable complex-valued function of uj, normalized so that 

cLu\f(uj)\ 2 = l, (2.5) 



where in this and the following integrals, the integration limits are — oo and oo unless other- 
wise specified; the role of negative frequencies will be discussed shortly. For the given mode, 
a^(oj) is the creation operator for excitation at an angular frequency uj. It is the adjoint of an 
annihilation operator a(uo), and the commutation relation between the two is 

[a(uj),a\uj')} = 8(u-uj'), (2.6) 
[a(uj),a(uj')] = 0. (2.7) 



This relation is a simplification appropriate to fiber modes of the commutation relation in [6]. 
The annihilation operator is required to satisfy the rule 

(Vw) a(u)\0) = 0, (2.8) 

which from the adjoint of Eq. (2.4) implies that 

a f \0) = 0. (2.9) 

From Eq. (2.6) follows the commutation rule for any functions / and g: 

[a g ,a*f] =JJ dudu' ' g*(u)f(u')[a(u),<J(u')] 
= J J dudJ g*{u)f{J)5{u-J) 

dug*(u)f(u), (2.10) 



where the asterisk denotes the complex conjugate. This can be written more compactly as 

[a g ,a\] = (gj), (2.11) 
where we define the inner product of functions 

(g,f)= [ dug*(u)f(u). (2.12) 



From this with normalized functions g = f, it follows that (/, /) = 1 by Eq. (2.5), and we 
see that and a/ are boson creation and annihilation operators satisfying [a/, a^] = 1. As a 
result, we can calculate the norm of a single-photon state a^\0) to be 

||4| )|| = ((Ola^lO)) 1 ^ 

(aja / + l)|0>] 1 / 2 = (0|0> 1 / 2 = l, (2.13) 



where the next-to-last equality follows from Eq. (2.9). For probabilities to make sense, we 
must require a finite inner product on quantum states, which requires normalizable states. This 
makes some bandwidth necessary: there can be no normalizable states at a single frequency. 
Because a^(uj)\0} has no norm, we call a\u) an improper operator. 



The hamiltonian operator for the mode (without zero-point energy) is 

H = h J du\uo\a ] {uj)a{uj). (2.14) 

Example 1: The expectation energy of a single-photon state a^\0) defined by the energy oper- 
ator of Eq. (2.14) is: 

(0\a f Ha\\0) = h(0\ J J J du x du 2 dco 3 f* (^ I )a(cu 1 )|cu 2 |a t (cu 2 )a(cu 2 )/(^ 3 )a t (^ 3 )|0) 

= h I du\to\ |/M| 2 , (2.15) 



where the second equality is obtained using Eq. (2.6) and integrating out the 5-functions. Thus 
if f(uo) is concentrated around some uo , one finds by this prescription a photon energy of 
about h\u G \. 

Commutation rules enable the calculation of probabilities of the form of Eq. (1.1). For 
instance, for a pure state, the right-hand side of Eq. (1.1) takes the form 

(0| J duj! - ■ ■ du n g(uj u . . . ^u^Poliaiiu^^liu^^biiu^^liuk), . . .)\0), (2.16) 

where Pol is a polynomial in creation and annihilation operators subject to ^-function com- 
mutation relations, and every annihilation operator acting from the left on the vacuum state 
|0) gives 0, as does every creation operator acting on the right of (0|. The standard method 
of evaluating such a probability is to use the commutation relations to put Pol in normal or- 
der. In the cases of interest, every power of every creation operator is paired with the same 
power of the corresponding annihilation operator, so that the only terms that contribute after 
the commutations have the form 

(0| J du 1 ---du n g(u 1 ,...,u n )Q(u 1 ,...,u n )\0) 

= J dw 1 ---dw n g(u 1 ,...,u n )Q(ui,...,u n ), (2.17) 

where Q is a sum of products of 5-functions, and the equality follows from Eq. (2.1). 
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A. Single-photon state spread over multiple modes 

Two modes a and b are called orthogonal if and only if [a, tf] = [a, b] = 0. Superpositions 
of single-photon states across orthogonal modes a and b have the form 

(c ia/ + 026^0), (2.18) 

where ^\ \cj\ 2 = 1. 

B. Single-mode, multi-photon states 

To construct the most general single-mode n-photon states, we introduce notation for multi- 
photon operators in which there are as many frequency variables as there are operator factors: 

h:a) n ^ J J < ^ UJl ' ' ' < ^ U ° n ^ UJl) ' ' ' ' UJn ) a ^ UJl ) ' ' ' a \ u n)- (2.19) 

Because the creation operators commute with one another, all that matters about h is the part 
of it symmetric under interchange of arguments, denoted 



S(u u ...,u n ) h(u u . . . , u n ) = — y~] h(u wl , u m ), 



vre5„ 



(2.20) 



where S n denotes the permutation group of order n. Thus we have 

h:a ]n = S(u 1 ,...,u n )h:a! n . (2.21) 

The normalization condition 

J' J du 1 ---duj n \S(uj 1 ,...,uJ n )h(uj 1 ,...,uJ n )\ 2 = l (2.22) 
assures unit norm for the single-mode n -photon state 

! (/i:a tn )|0>. (2.23) 



'71! 

The adjoint works according to 

[(/i:a t ")|0)] t = (0\(h*:a n ). (2.24) 



1. Inner product of two single-mode, multi-photon states 



Consider two modes that need not be orthogonal, such as two linearly polarized modes a and 
b with an angle 9 between them, so the commutation relation between them is [b(u), al(u')] = 
cos 9 5(uo — uo'). For symmetric functions h of n arguments and h! of m arguments, the usual 
manipulations show that 

-l=(0\(h'*:b m )(h:a^)\0) = (S^, . . . , u m ) h', S(u u . . . , u n ) h) cos n 9, (2.25) 
vn!m! 

where we define the inner product of the multi- variable functions h and h! as 

{0, if h and h' are unlike in number of arguments, 
(2.26) 
/• • •/ dwi--- du n h'n(ui, . . . , u n )h n (ui, ...,u n ), otherwise. 

Thus for h — h' and h normalized per Eq. (2.22), the inner product of the two n-photon states 
is just cos" 9. 

2. Energy of single -mode, n-photon state 

For h normalized per Eq. (2.22) and the energy operator H of Eq. (2.14), the expectation 
energy for the state n\~ l l 2 (h:a vn )\fy is 

^-(0\(h*:a n )H(h:a^ n )\0) = nhu h , (2.27) 
n\ 

where we define an average angular frequency 

^ft= f J "J dwi - ••dwn\ui\ \S(u 1 ,...,u n )h(u; 1 ,...,u n )\ 2 . (2.28) 

C. Broad-band coherent states 

As a special case of a two-photon state, one can choose any normalized function / and 
define h(ui,U2) = 2~ 1 / 2 f(u 1 )f(u2) to produce a two-photon state 2 _1 / 2 (aj) 2 |0). This imme- 
diately generalizes to higher powers of a^, leading to a 'broad-band' coherent state. From 
Eq. (2.10) we see that ay and aj are boson creation and annihilation operators satisfying 
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[a/,ajf] = 1, so that Louisell's discussion of coherent states [7, Sec. 3.2] applies to them. 
Hence we can define n-photon /-states by powers of acting on the vacuum state n times, 
producing the normalized state (n!)~ 1 / 2 (aj.) n |0). Coherent /-states can be defined as super- 
positions of these states, just as in Louisell [7, p. 104]: 

\a,a f ) d = f e W (~\aAjr(n\)-\aa\) n \0) 

^ ' n=0 

= exp ^-i |a| 2 ^ exp(aa\)\0). (2.29) 

This is an example of a calculation that proceeds exactly as if frequency dependence were 
collapsed, so that aj works like a simple oscillator annihilation operator. In contrast, frequency 
dependence matters in the commutation rule 



<kjg*(u)f{u), (2.30) 



which follows from Eqs. (2.6) and (2.7). 

To calculate the energy of the coherent state, note first that the commutation rules Eqs. (2.6) 
and (2.7) imply for the coherent state defined by Eq. (2.29) 

a(cu)\a,af) — af(cu)\a,af). (2.31) 
From this and Eqs. (2.14) and (2.5) follows the expectation energy 

(a, af\H\a, a/) = ftwf\a\ 2 , (2.32) 
where we have defined a mean absolute angular frequency 

u f = J du\u\\f{u)\ 2 . (2.33) 

D. General state 

The general state is a sum of terms, not necessarily normalized, each of the form 

j 

i^)=/:n a ri°)- < 2 - 34 ) 



Symmetries under interchange of variables are best expressed in a notation intermediate be- 
tween that of writing out the integrals and the compact "colon" notation. We expand the short- 
hand / : af n to /(u>) : a) n (<jj), understanding that with an exponent n involved, (u>) is short for a 
list of n frequency variables and a) n (u>) is short for af(ui) ■ ■ ■ at(u n ). Written in this notation, 
Eq. (2.34) becomes 

j 

|V> = /(wi, , wj) : J] ( Wj )|0); (2.35) 
i=i 

where (ut,) is a list of n,- frequency variables. Symmetry under interchange of creation opera- 
tors implies, as in Eq. (2.21), 

/(w):a^(w) = [5(w)/(w)]:o^(u;), (2.36) 

where the symmetry operator S is defined in Eq. (2.20); this generalizes to 

j j 
/(«!, , wj) : aj n ' ( Wi )|0> = • • • . . . , uij)] : JJ a*' («,-) |0). (2.37) 

j"=i j"=i 

Sometimes we write such expressions a little more compactly, using the convention that uj n = 
(u>i,...,wj) as 

j 



/K) : ns K K)i°) 



L \i=i 



j'=i 



We can distinguish one set of modes from another by replacing some of the aj by other labels, 
such as bj. For example, a-modes can refer to Alice while 6-modes refer to Bob. This leads to 
a general term in a state expansion of the form 



J \ / K 

J n 3 



/K,a>„): ( nV'K) n^^) ) l°>- ( 239 ) 
vj=i / \fc=l 



E. Density matrices and traces 

A density matrix is a sum 



p = ^w n \*(n))(*(n)\, (2.40) 
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where the \^(n)} are unit vectors, w n > 0, and J2 n w„ = 1; or an integral 

p = Jw(u)du\V(u))(V(u)\, (2.41) 

with the \ty(u)) unit vectors, w(u) > and J duw(u) = 1. More generally, a density matrix 
can be any convex sum of these discrete and continuous types. 

We need the trace of a density operator multiplied by a bounded operator M. Although 
in the context of infinite-dimensional spaces the trace is sometimes defined only for positive 
operators [8,9], we want to apply a trace to terms that occur when a pure- state density operator 
is expanded (i.e., |*)(*| = Y. m , n c m c* n \^ m )(^ n \. For this we define 

Tr[(|^>(0|)M] = (0|M|^>, (2.42) 

and define more general traces by linearity. We then find 

Tr(pM) d = J^iu n (*(n)|M|#(ra)) (discrete case), 

n 

Tr(pM) d = Jduw(u)(m(u)\M\m(u)) (continuous case). (2.43) 
F. Partial traces of light states 

Probabilities for detection of light states often involve partial traces, defined for an operator 
on a tensor-product space Ha^'Hb by linearity from the following special case. For |o), \a!) G 
Ha and \b') eH B , 

Tr B (|a)|&><6'|(a'|) = (W»|a)<a'|, 

Tr A (\a)\b)(b'\(a'\) = ((a'\a))\b)(b'\. (2.44) 

From this it follows that for any bounded operator Ma that acts only on Ha, we have 

Tr A [(\a)\b){V\(a'\)M A ] = ((a'\M A \a))\b)(b'\. (2.45) 

Here is an example. If a detection operator is of the form M A ,j ® 1b, then 

T*AB[pAB{M A ,j ® l B )} = Tv a (pMa,j), (2.46) 
13 



p' = Tr B (p AB ). (2.47) 

The partial trace over 6-modes of any light state p is calculated from the commutation relations 
applied to the inner product of 6-mode factors in the usual way, supported by the notation 
developed in the preceding subsection. For example, one can deal with the case of a single 
a-mode and a single 6-mode as follows. Suppose we have 

\i/> mn ) = jdudug(u,u)a^{u)b^{u)\0), (2.48) 

where we abbreviate u = u m , uj = a>i, . . . , u> n , a^ m (oj) = a\uji) ■ ■ ■ a)(uo m ), etc. 



Similarly, suppose 

\<j> m , n ,) = Jdu>'duj'h(u>\uj')a^ m \u>'p n '(Cj')\0). 
Then, for an operator Mj(a) that commutes with b(u>) and b\u)), 

Tr ab [(|^ m „)(0 m '„'|)M j (a)] = TrJp'M^a)], 

where 

p = Tr 6 (|^mn)(0m'n'|) 



(2.49) 



(2.50) 



J dujduj' J du>du>'g(u>, u)h*(ut', u) (0 b \b n ' (a>> tn (u>) |0 6 ) 



a tm M|0 a )(0 a |a m V) 



= n\5 n „> Idujdu}' 



a^(u)\0)(0\a m '(u'), 



(2.51) 



' duj g(uj, uj)S(uj)h*(uj' , uj) 

where the second equality follows from Lemma (B16) of Appendix B applied to the 6-modes. 
Similarly, one calculates for the trace over the a-modes 

Tr a (|Vw)(0m'n'|) 



dujduj' 



dujduj'g(uj,uj)h*(uj' ,Cj)(O a \a m ' (uj')a^ m (uj)\O a ) 



tf n (u>)\0 b )(0 b \b n '(u>') 



m\5 mrn i Jdujduj' Jdujg(uJ,uj)S(uj)h*(uj,uj') 



6 tn (w)|0)(0|6 n '(w')- 



(2.52) 



The more general case involves complications, but here it is. Let n = J2j=i n j an d m = 
J2k=i m k an d similarly for primed quantities; let 



\j=l / \k=l 

J' \ / K' 

K 



l^> = ^'K„<,): ( II VK) II 6 * ) l°>- 

Kj=l J \k=l 



(2.53) 
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On expanding the 'colon' notation and abbreviating using duj n for du>idu> 2 ■ ■ ■ dujj and duj r 
for duj\duji ■ ■ ■ duj K , one finds 



= / doj m du>' m , 



j' 



J du) n dv' nl h(u> n , 6) m )h!*(u)' n ,, Gi'J) 



n^(«*)l io fo ><o fo j fnc'n^)) • (2.54) 

^fc=l / \fc=l / 

[The last equation follows from Lemma (B15) of Appendix B.] Because the 6-creation opera- 
tors commute among themselves, as do the b- annihilation operators, there is one more symme- 
try: 



= <5 j jf ^^nj\8 njn ^j JdCb m duj' m , j Jdu:, 



K 



iis(u; k )\ n%) ^ 

n ■ 



\k=l 



\k=l 



\j=l 



K 



K' 



n4 mfc (^))io b )(o 6 i(ncn^: 

(2.55) 



\k=i 



\k=i 



G. Bi-photons: excitation in each of two orthogonal modes 

For a state that exhibits a single photon in each of two orthogonal modes, a and b, whether 
in a single fiber or in different fibers, the general form is 



(h:aW)\0), 
where we extend our notation by defining 

(h:a^tf) = f / / du du h(u, u))a^(u)b^(u>), 



(2.56) 



(2.57) 
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with the normalization requirement that 



1 



// 



du duj \ h(u,cu)\ 2 . 



(2.58) 



Note: the normalization requirement rules out h(u, Cu) of the form f{uj)5{uj — uj). 

Other states, including multi-mode, multi -photon states will be introduced after we have the 
machinery of projections. 

3. PROJECTIONS 

To deal efficiently with polarization-entangled states, we need to characterize subspaces of 
states by the projections that leave them invariant. These are provided here. It is instructive to 
compare and contrast the projections for the function space of states used here with the case 
of a single oscillator with its basis states \n), n — 0, 1, . . . [7]. With our range of frequencies 
there are countless 1-photon states (and countless n-photon states), as described in Sec. 2B, 
in contrast to the single n-photon state \n) for given n of an oscillator; however, there is a 
one-to-one correspondence between a set of projection operators for photon number and the 
projections \n)(n\ for the oscillator. 

For both the single oscillator and the states that are characterized by functions of frequency, 
one writes |0)(0| for the vacuum projector. To the oscillator projection |1)(1| corresponds the 
projection that leaves invariant single-photon states while killing all states of more or fewer 
photons: 



It is easy to check that for any / one gets Pi (a) (a^ | 0)) = a | 1 0), and that for any of the n-photon 
states \ip n ) with 1, one finds Pi (a) \ip n ) = 0. Similarly the projection for two-photon states 
is 




(3.1) 




(3.2) 



The general case is 



cLui ■ ■ ■ du) n a}{uji) ■ • • a) (cu n ) |0) (0|a(c^i) • • • a(cu n ). 



(3.3) 
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It follows that any state in the space of superpositions over all n is unchanged by the operator 
obtained by summing over all P n (a), and hence this sum is the unit operator 



X>» = 1- (3.4) 

n=0 

Equation (3.3) implies a relation among the projections for adjacent values of n 

P n +i(a) = — -J— j- J duja ] (u)P n {a)a(u)). (3.5) 

Although this mirrors the situation for a single oscillator, the projections P n (a) cannot be 
expressed as the outer product of a vector with its adjoint; indeed the n-photon subspaces are 
infinite dimensional. 

From the commutation Eqs. (2.6) and (2.7) it follows that 

m m 

P n (a) n« f K)|0) = 6^11^(^)10); (3.6) 

3=1 3=1 

furthermore, this holds if one inserts into both sides of the product the same polynomial in 
operators for modes orthogonal to a. Also from Eq. (3.6) follows the relation 

P n (a)(h:a^ m )\0) = <5 n , m (/*:a tm )|0>. (3.7) 



A. Action of single-mode projections on multi-mode states 

In the context of two orthogonal modes a and b, the symbol P n (a) is re-used as shorthand for 
P n (a) <g> lb, where l b is the unit operator on the 6-factor of a tensor product. Two projections for 
differing values of n for a given mode are mutually orthogonal, but not when one projection is 
for one mode and the other projection for another; e.g. P 1 (a 1 )P s (a 1 ) = but Pi(a 1 )P 3 (a 2 ) ^ 
0. Similar remarks apply to the case of more than two mutually orthogonal modes. 



B. Multi-mode n-photon projector 



To study polarization we will need the idea of an n-photon state that can be distributed over 
two modes. We denote the projector for n photons distributed arbitrarily among modes a x and 
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a 2 by 

n 

P n (a 1 ,a 2 ) 

fc=0 

Note: All these projectors can be subdivided into "+" and "— " parts, e.g. 

P n (a) = P n , + (a)+P n ,4a). (3.9) 

C. Number operator 

Similarly to the single oscillator case, one constructs an operator that has as its expectation 
value for a given state the mean photon number for that state. This is the number operator for 
the mode under discussion: 

oo 

N(a) d ^J2 nP n( a )- ( 3 - 10 ) 

n=0 

By use of Eqs. (3.5) and (3.4), this is transformed into the more convenient form: 

N(a) = } nP n (a) = y^n— / duo a) (uo)P n _i(a) a(u) 
i i n J 

n=l n=l 

= J dua\u) ^^JP„_i(a)j a(u) = j duo a\u)a(u) . (3.11) 
It is also useful to express the number operator for the "+" and "— " directions: 

/'OO 

N±(a)= duoa ] ± {uo)a±{uo). (3.12) 
Jo 

For a space (or subspace) spanned by two modes a x and a 2 , define 

OO 

N(a 1 ,a 2 ) = y^ j nP n (ai,a 2 ). (3.13) 

?1=0 

As sketched in Lemma (B21) of Appendix B, a calculation similar to that for one mode shows 

N(ai,a 2 ) — J du[a\(u)ai(uj) + a 2 (uj)a 2 (uj)]. (3.14) 

Here again, we can pick out directions, as in Eq. (3.12). There is no difficulty in extending to 
more than two modes. 
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4. LOSS AND FREQUENCY DISPERSION 



We partition the modes to be analyzed into "desired modes" and "extraneous modes," as 
illustrated in Fig. 2. By considering coupling of desired modes to extraneous loss modes, one 
can model a wide variety of loss mechanisms. The choice of model for loss is tied to the choice 
of model of detection. 

To model loss of a mode to be detected by a binary detector, neglecting memory effects in 
the detector, for now we succumb to the charm of simplicity in the approach of Mandel [10, 
p. 640] . In this approach a mode b' prior to loss is related to a mode b after loss by the following 
equation for the respective annihilation operators: 



where c expresses an undetected mode into which all the fe'-energy spills except for a fraction. 
In the case where ?7 loss is independent of frequency, this un-lost fraction is just |r] loss | 2 . 

An alternative is to consider the coupling of desired modes to a heat bath. This leads to 
so-called master equations. The simplest form that expresses the essential features makes the 
Markhoff approximation [7, p. 347, Eq. (6.2.61)]. With further simplifications, including that 
of a zero-temperature heat bath, we obtain the time behavior of a density operator that at time 
t expresses the state at 1 0) to be, in the Schrodinger picture (SP), a density operator of the form 



[Need to work this out for multi-photon states.] 
A. Loss cannot evade "no cloning" 

I believe (and need to check) that coupling of desired modes to thermal modes cannot make 
two states more distinguishable. More formally, suppose the desired modes are A and the 
extraneous modes are B, and the in-state is expressed by either some tensor product pA,m,i <8> 



b'(u) = r] loss (uj)b(uj) + (1 - |^7ioss(^)| 2 ) 1/2 c(^), 



(4.1) 




(4.2) 
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PB,m or by pA,in,2 <S> PB,in- After a unitary evolution U acting on the tensor-product space of A 
and 5, the out- state is U PA,m,j <S> PbU\ with j = 1 or 2. The reduced density matrix for A is 
obtained by the partial trace over B: 

pA,out,j = Tr B (Up A ,m,j ® PbU ] ). (4.3) 

(Note that, unlike traces, partial traces of a product depend on the order of the factors.) My 
guess is that in all cases 

Tr(p^ util pY 2 ou t, 2 ) > MpIIAIJ- (4-4) 

For this it is necessary and sufficient to prove for all density operators p\ and p 2 acting on 
the tensor-product space of A and B that 

Tx A [(Jx B Pi) 1/2 (Tr B p 2 ) 1/2 ] > Tr(pfpf ). (4.5) 



5. LOCAL QUANTUM FIELDS 

Improper operators introduced so far have been integrated over frequency to produce proper 
operators. We want also to construct proper operators by integrating over space and/or time, 
which leads us to take Fourier transforms, as outlined in Appendix D. In analogy with quantum 
electrodynamics for propagation in vacuum, we introduce (improper) local annihilation field 
operators in the Heisenberg picture: 

1 r°° 

a+(x,t) = -= / du;a + (u;)e- i[u}t - k ^ x] , 

a-(x,t) = -= \ duo a(uj)e~^ t+k ^ )x \ (5.1) 
\/2Wo 

where k{uo) is an experimentally determined propagation factor, defined so that sgnk(uj) = 
sgncu. We assume 

1. fc(-w) = -jfe^). (5.2) 

2. lim^^oo > o. (5.3) 
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3. *P- > 0. (5.4) 

4. Kramers-Kronig relations connect this k(u) to the loss coefficient ^(uo) of Eq. (4.2) 
[11]. 

The operator a±(x, t) is 'improper' in that it takes a normalized state to an unnormalizable 
state; a proper operator can be defined by averaging a±(x, t) over a spacetime region. (This av- 
eraging takes place automatically in time-dependent perturbation theory [7, p. 257].) Drawing 
on Eq. (2.2), we can define 

a(x,t) = ■ a + (x,t) +a-(x,t) = -= / dcu a(w)e- i[M *- fc(w)x] , (5.5) 

V27T J-oo 

and from Eqs. (5.1) and (2.6) follows the commutation relation 

[a(x, t), at(u)] = exp{— i[\u\t — k(u)x]}. (5.6) 
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A. Temporally local hermitian fields 

From the non-hermitian operator a(x, t) can be constructed the two non-commuting hermi- 
tian quadrature operators 

q(x,t) = f a)(x, t) + a(x, t) 



duo {a t He i[|w|t - fc(w)x] + a(tu)e- l[Mt ~ k ^ x] }, (5.7) 



oc 



\f2n J 

p(x,t) = f i[a\x,t) — a(x,t)] 

= -== I dw{a t (w)e i[|w|t - fc(w)x] -a(cu)e- i[l " l *- fe( " ):El }. (5.8) 

Note t) is like Louisell's voltage operator. Fourier-transforms show 

/oo />oo 
dt a±(x, t)a±(x, t) — / duo a±(u)a±{u). (5.9) 
-oo JO 

This will be used in constructing operators to model detection. 
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B. Time, space, and dispersion 



Because the variation of k{u) with u is non-linear, the operator that is convenient for space 
localization must differ from that which is convenient for time localization. Consider op- 
erators such as the hamiltonian and the number operator (shortly to be introduced) of the 
form O = f du f(u)a\u)a(u). In the case of the number operator defined as N±(a) = 
J °° du a}±(u)a±(u) = J dt a±(x, t)a±(x, t), the local operator a±(x, t)a±(x, t) acts as a kind 
of density in time for N±(a); however, the space integral of a±(x, t)a±(x, t) is something else: 

/OO /"OO poo 

dxai(x,t)a±(x,t) = / dw x \ du 2 a^u^^y^'^d^u^-k^)], (5.10) 
-oo J0 JO 



which, with the relation 



S[k(u) - k(w')] 



dk{uj) 



du 



5{uj — u'), 



becomes 



[uo 



dk{uj) 



du 



/oo poo 
dx a±(x,t)a±(x,t) — / du a±(uj)a±(- 
-oo Jo 

Note that the "extra factor" l/\dk/du \ expresses a group velocity [12]. 
With this in mind, we construct an energy density operator a'±(x, t) by 



a±(x,t) 



def 



du ( U 



dk{u) 
du 



1/2 



a±(u)e 



One can check that this has the property 

/OO pOQ 
dx a±(x,t)a'±(x,t) — / du \u\a±(u)a±(u) — H/h. 
-oo Jo 

There are lots of other possibilities. For instance, define a field 

def 1 



a(x, t) 



then 



H/h 



du 



. _ + . . da(x,t) 
dxa'(x,t) — , 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



so that the two distinct fields a and a both enter. 
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In order to deal with finite fibers, we explore the operator obtained by making the limits of 
integration over x finite; however, the resulting integral is no longer independent of t. Further, 
while the integral obtained from putting finite rr-limits in (5.14) is at least hermitian, even this 
modest property fails for the integral obtained by putting finite rr-limits in (5.16). Nonetheless, 
in some cases such a truncated operator can be a useful approximation to the hamiltonian. 

C. Projections in terms of local operators 

It will be interesting to examine approximations to projections in terms of operators that 
are, so to speak, confined in time and space. For reference, here we express some projections 
in terms of local operators. From Eq. (3.1) and the inverse Fourier transform in time of Eq. 
(5.1), one computes 

/oo 
dta^(x,t)\0){0\a(x,t) (5.17) 
-oo 

(independent of x). Similarly, from Eq. (3.2) one computes 

/oo ;>oo 
dt / dt'a}(x,t)a^(x',t')\0)(0\a(x,t)a(x',t'), (5-18) 
-oo J — oo 

and one can keep on going. This procedure also works to express the number operator of Eq. 
(3.1 1) in terms of local operators: 

N(a) = f dwa\u)a{u)= f dta\x,t)a{x,t). (5.19) 

J J — oo 

6. SCATTERING MATRIX 

Couplers and other networks can be analyzed in terms of a scattering matrix. For this, one 
supposes that there is time t\ before which the network acts as a set of uncoupled fibers and 
perhaps a quantum memory uncoupled to these fibers; one supposes a later time t 2 after which 
the interaction is over, so that again one has a set of uncoupled fibers and a quantum memory 
in isolation. Prior to t±, the light is modeled by an in-state as an integral over frequency of 
a polynomial in in-mode creation operators acting on the vacuum. For time t > t 2 the light 
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is modeled by an out-state, again as an integral over frequency of a polynomial in out-mode 
creation operators acting on the vacuum. 

For example, to analyze a coupler, we typically express state preparation in terms of in-state 
creation operators while we express detection in terms of out-state creation operators. Hence 
calculating the probability of an outcome calls for expressing in-state operators in terms of 
out- state operators, or vice- versa, according to whichever is more convenient. 

The use of scattering theory in this context is to relate the creation operators for in-state 
modes to creation operators for out-state modes, and this relation is constrained to be a unitary 
transformation. Because of unitarity, commutation relations are preserved and so are inner 
products. A complication is that unitarity holds only for a transform over all the modes in- 
volved, "extraneous" as well as "desired." 

In general the creation operator for an out-mode at one frequency depends on creation op- 
erators for in-modes at all frequencies; i.e., the operator mixes frequencies, as in parametric 
down conversion. For linear networks we get the major simplification that out-mode operators 
for a given frequency uo are unitary transforms of in-mode operators for the same frequency. 

After a scattering transformation is obtained, either by guessing or by calculating it from 
some more detailed model, one can check whether the in-states and out-states are consistent 
with particular values of t\ and t 2 . (For a network with strong internal reflections, the interval 
t 2 — ti must be large enough to allow reverberations to die out.) 



A. Network without frequency mixing 

For the special case of a network that is linear in that it mixes no frequencies, if the extrane- 
ous modes are accounted for, then the out-state annihilation operators aj jOUt (uj), j — 1,2, ... , 
are some unitary transform of the in-state annihilation operators. That is, let 



^out 



Ol,out 
°2,out 



(6.1) 
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and similarly define a in . Then there is some frequency-dependent, unitary matrix U (uS), having 
dimension equal to the number of out-modes (by assumption equal to the number of in-modes), 
such that 

aout(^) = U(u)a in (u). (6.2) 

Taking hermitian conjugates and multiplying the result by U (uo) on the right gives a relation 
between the row vectors a} ont and a\ n : 

oL = 4^H- (6-3) 

To understand the application of scattering formalism to modeling fiber interactions illus- 
trated in Fig. 1 , first model light pulses as quantum states expressed by creation operators acting 
on vacuum. We think of an interaction that, to some approximation, starts at t\ or later and is 
complete by t 2 or earlier. Prior to ti, the incoming light is modeled as involving only in-state 
creation operators. After t 2 , the light is modeled as a different state, involving only out-state 
creation operators, related to the in-state creation operators by an equation of the form (6.3). 

7. POLARIZED AND ENTANGLED LIGHT STATES 

It is possible to make a fiber, for instance with an elliptical cross section, that propagates 
only a single mode; however, the fibers used in the Quantum Network all carry two polarization 
modes. We follow convention by misnaming fibers, so that when we say a 'single-mode' 
fiber we mean a fiber that propagates not a single spatial mode a but one that propagates two 
polarized modes a\ and a 2 , mutually orthogonal in that the [ai,c4] = 0, with propagation 
constants k a i(uo) and k a2 (uo), respectively. This means that for a network (as in Fig. 1) with n 
'single-mode' fibers there are 2n in-modes and 2n out-modes, so that, neglecting extraneous 
modes, the single-frequency scattering matrix U (lu) has dimension 2n. 

(For special materials, one may need to separate out the two directions to allow k aj+ ^ 
k a j~.) Note that a linear superposition of the two modes with different propagation constants 
has no definable propagation constant. 

In describing modeling approaches for fibers that support two polarizations we make (and 
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* (a>) x = k b .((a) 

FIG. 3: Splice in fibers of different type. 

try to state) various simplifying assumptions; throughout we assume frequency conservation. 
Taking polarization into account doubles the number of modes relative to the number of so- 
called single-mode fibers, for the simple reason that these are misnamed. 

For a polarized fiber having orthogonal modes a% and a 2 , the form of a single-photon state 
is [see Eq. (2.18)]: 

1 1 -photon) = (cidij + c 2 a 2)g y\0), (7.1) 

for any normalized / and g, assuming numerical constants \ ci\ 2 + |c 2 | 2 = 1. The same operator 
generates the coherent state of the form of Eq. (2.29): 

|a,(ciai,/ + C 2 a 2i9 )) = e -|°l 2 /2 e a(c 1 a ll/+ c 2 a 2 , B )t| 0) 

= \c x a, a x j) <g> \c 2 a, a 2 , g ), (7.2) 

where the second equation follows from exp[(ci<2ij + c 2 a 2)9 )t] = exp(ciaj j) exp(c 2 a 2j9 ) 1 ", 
which is a consequence of the commutativity [a\, a 2 ] =0. 

[**The whole business seems to assume a reference position x = 0; other values of x 
involve x -dependent phase factors.] 

A. Fiber splice (without extraneous modes) 

Consider the situation shown in Fig. 3, where a fiber a on the left is spliced to a fiber b 
on the right. We suppose the splice is centered at x — 0. For such a heterogeneous spliced 
fiber, the local field operator a(x,t) cannot satisfy a field equation that exhibits translational 
symmetry; instead, the field equation stems from a hamiltonian that has a change at the splice; 
the situation is reminiscent of a potential problem in quantum mechanics. 
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Without trying to model the details, and neglecting coupling to extraneous modes, the effect 
of the splice is to convert in-mode operators to out-mode operators, as expressed by a unitary 
matrix U (u). Recognizing polarizations, we have modes ai and a 2 in the a-fiber and modes b\ 
and b 2 in the b- fiber, with propagation constants k ai , k a2 , fc bl , and kb 2 , respectively. This makes 
the four in-modes a j+ and bj-, j = 1, 2, along with the four out-modes a 3 -_ and bj + ; hence if 
we neglect loss, the scattering matrix U (lu) is 4-by-4: 

[a[_(u), b[ + (u), 4_H, = [&I-M, <A + (u), b\_{u), a\ + {u)]U\u). (7.3) 

Under the simplifying assumption that the splice has no coupling between polarizations, the 
equation factors into two 2-by-2 pieces, one piece for each polarization. With U (uo) defined this 
way, the perfectly homogeneous situation in which k aj (uj) = k bj (u), thus making the splice 
invisible, corresponds to U(uo) = 1. 

(To account for loss, we can add extraneous dimensions and then trace them out, thereby 
getting a matrix that is not unitary, so that the out-power can be less than the in-power.) [*For 
single-frequency matrices, we use power rather than energy, because energy is definable only 
for non-zero bandwidth.] 

B. Coupler 

Fiber couplers, analogous to beam splitters, have four fibers, aj,bj,Cj,dj, with j = 1,2 
(see Fig. 4). Each of the eight modes comes in two directions, "+" and "— ". Without loss, 
U(u) G SU(8). 

C. Entangled states 

Whether entangled in frequency or in polarization or in both, entangled states have to do 
with tensor products of vector spaces. Let a vector space V be a tensor product of vector 
spaces V a , where a ranges over some index set. Relative to this factorization, a generic vector 
in v E V is a sum of tensor products of vectors v a G V a ; unless it can be written as a single 
product, not a sum of products, v is called entangled (relative to the factorization). 
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FIG. 4: Fiber coupler. 



We stress relative to a factorization because vector spaces for light states can be factored in 
more than one way, and to speak sensibly of 'entangled states' one must know which factor- 
ization is meant. The vector space for a single mode involves tensor products over subspaces 
for frequency bands; this factorization is relevant to frequency-entangled states. Within any 
frequency band, there is an infinite tensor product over single frequencies, which will not be 
used here in speaking of 'entanglement'. Polarization-entanglement involves a vector space 
for several modes, factored by mode. 

Frequency-preserving, unitary transformations of the single-frequency creation operators 
carry an in-state that is a tensor product of broad-band coherent states to an out-state that is 
also a tensor product of broad-band coherent states. This is a special property; indeed the only 
case I know in which it works is that of coherent states. [* Work out proof] (In generic cases, 
a tensor product of single-photon in-states is carried to an entangled, multi-photon state.) This 
makes broad-band coherent states of special interest as an (overcomplete) basis for studying 
light in the context of frequency-preserving splices and couplers. See Appendix E. 

D. Polarization-entangled states 

A polarization-entangled state involves four modes; we think of two orthogonal modes a x 
and a 2 in a left fiber to Alice and two orthogonal modes b\ and b 2 in a right fiber to Bob. (One 
can think of 1 as horizontal and 2 as vertical.) Then there are four creation operators a\(u>), 
a\{uj), b\(uj), and b\{uj)\ everything commutes except that [cij(u;), aj(u/)] = S it j5(uj — u/), 
where 5 i: j = 1 if j = i and otherwise is 0. 
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The defining property of a state \(f>) characterized by a single a-photon combined with a 
single 6-photon is 

P 1 (a 1 ,a 2 )P 1 (b 1 ,b 2 )\<f)) = \<f>). (7.4) 
A necessary and sufficient condition for this is that |0) have the form 

2 

I0> = E^m(M&J)|O>, (7.5) 

where 

2 

where we use the definition stated in Eq. (2.57), and the four otherwise arbitrary functions h i: j 
satisfy the normalization condition (in which no symmetry is assumed): 

~.\|2 



du dw \hij(u), d))| =1. (7.7) 

To think about quantum key distribution, we will need to consider states of both fewer and 
more photons. A state with m photons in a-modes, j of which are in mode a±, and with n 
photons in b- modes, k of which are in mode bi, has the form 



a{ a r%br k m 



y/j\(m — j)\k\(n — k)\ 

= — f - f dwidwi • • • duj n dw n h(ui, u>i, . . . , u n , u n ) 

y/j\(m-j)\k\(n-k)\J J V ^ 

x aj(wi) • • • a{(uj)oi(u j+1 ) ■ • ■ 4(w ro )6}(a>i) • • • b\(Cb k )b\(Cb k+1 ) ■ ■ ■b\{Cu n )\V>) . (7.8) 

There is no loss of generality in requiring h to be symmetric under each interchange ujj <-> uj k 
for which ujj and uo k pertain to the same mode. When h with this symmetry is normalized as in 
Eq. (2.22), the state defined by Eq. (7.8) has unit norm. The most general state is a weighted 
sum of such states. 



8. DETECTION 



A light detector as used in an experimental setup will be modeled by a positive operator- 
valued measure (POVM), {Mj}, so the probability of outcome j is Tr(M,p), where p is a 
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density operator, and J2j Mj = 1 with Mj > 0. In the simplest models to be discussed here, 
the density operator is for the light; in more complex models discussed elsewhere [2], the 
density operator can be for a pure or mixed state of not only light to be detected but also probe 
particles. One could also include the generation of backward-propagating light, but here we do 
not get into this level of complication. 

A. Simple examples 

Example 2: As a first theoretical example, a detector as narrow-band as possible that will regis- 
ter outcome 1 with certainty, given the one-photon state a* |0), has M 1 = a g \0)(0\a g . This 'fil- 
tered' detector discriminates well against any other one-photon state aj-|0) if / is unlike g; i.e., 
the probability for outcome 1 for a state a\\0) is Tr(aj|0)(0|a 9 a}.|0)(0|a / ) = | (0|a 3 aj.|0) | 2 = 
\jdug*{u)f{u)\\ 

Given a set of mutually orthogonal functions gj, j = 1, . . . , N, there can be a detector with 
N possible outcomes, with Mj = a g .\0)(0\a gj . Thus a single-photon detector is certainly not 
restricted to giving a yes-no outcome, and discrimination among different single-photon states 
is possible. Herein lies a caution to designers who hope that distinct lasers generate distinct 
states that differ in polarization. 

Example 3: An example of a model of detectors preceded by filters used for joint detec- 
tion of a two-mode state is Mi = a^b g \0){0\b g af. Then outcome 1 would be described as 
'finding the state Of\0) in measuring mode a jointly with finding the state b g \0) in measuring 
mode b.' For the two-mode state (2.56), the probability of this outcome is readily calculated to 

be \JJ dudu' f*(u)g*(u')h(u,u')\ 2 . 

Example 4: An interesting and perhaps novel application of probe particles involves detection 
at two coordinated locations using probe particles that have previously become entangled. This 
can produce an example that swaps the detection operator M 1 and the state of Example 2, leav- 
ing the probability invariant. I.e., the same probabilities and outcomes arise from measuring 
an entangled state with an unentangled detector as from measuring an unentangled state with 
an entangled detector [2]. 
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B. Model of APD detector for quantum cryptography 



For multi-photon states, a simple approach is to model a detector as a POVM for a state 
involving only the light to be detected. This assumes that memory effects in the detector, 
such as those often attributed to trapped carriers in photo-diodes, are insignificant (for instance 
because the network design enforces enough hardware dead time); it also assumes teetering 
in the detector plays no significant role [2]. For simplicity, here we add the assumption that 
the detector and its pulse-shaping circuitry choose an outcome of 'yes' or 'no' without any 
additional detail. With these assumptions, all one has left is dark-count and efficiency; one 
obtains a class of detector models that respond to n photon states according to 

Pr (Detect) n-photon state) = d n , (8.1) 

with < d n < 1. With the additional assumption that the detector is flat in its frequency 
response over the bandwidth of the light pulses to be detected, the detection operator that 
generates these probabilities is 

oo 

M 1 (a) =J2 d n Pn,+(a), (8-2) 

n=0 

where P n (a) is defined in Eq. (3.3) and we use a + (uo) in place of a(uo), so the integrals are over 
positive frequencies only. [For a detector of both polarization modes a x and a 2 of a fiber, one 
has instead, 

oo 

Mi(ai, a 2 ) = ^ rf n^n,+ (oi, a 2 ), (8.3) 

n=0 

with P n>+ (ai, a 2 ) defined by Eq. (3.8) with a + (oo) in place of a(u).] 

To model the use of avalanche photo-diode (APD) detectors in the DARPA Quantum Net- 
work, we explore in detail a specialization of this type of model. We will use this model in 
Sec. 9 to study the variation in detection statistics as the energy of transmitted pulses is raised 
above that of a single-photon state. 

An APD detector for an optical fiber responds to both polarization modes of the fiber; how- 
ever, by placing a polarizing beam splitter before the detector, one can effectively eliminate 
the light from one mode. We assume this case, so that only one polarization mode is relevant. 
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Then we suppose that the probability of the detector responding to any single-photon state is 
rjdet, assumed independent of frequencies over the range of frequencies relevant to the light 
state. We need to model the probability of detecting multi-photon states. By the assumption 
already made, the probability of failing to detect any single-photon state is 1 — r] dct . Assume 
that the probability of failing to respond to an n-photon state is just the single-photon failure 
probability raised to the n-th power: (1 — r] dct ) n . Then the probability assigned by this model 
to registering a detection, given an n-photon state and temporarily ignoring dark counts, is 

Pr (Detect) n-photon state) = 1 - (1 - r] dct ) n . (8.4) 

To allow for dark counts we replace this by 

Pr (Detect | n-photon state) = 1 - (1 - p daik )(l - r] dct ) n . (8.5) 

To deal with states that are superpositions over varying numbers of photons, we recall that the 
operator for an ideal detector that responds always to any n-photon state, with no probability of 
detection for a state that has no n-photon component, is the projection P n> ± defined in Sec. 5. 
The following detection operator Mi (a) invokes P n> ± to provide the desired probabilities for a 
detector of a single a-mode [13]: 

oo 

M x (a) =^[l-(l-p dark (a))(l-r ?dct ) n ]P n , + (a), (8.6) 

n=0 

where P n (a) is defined in Eq. (3.3) and we use a + (uo) in place of a(uo), so the integrals are over 
positive frequencies only. [For a detector of both polarization modes a x and a 2 of a fiber, one 
has instead, 

oo 

Mi(ai, a 2 ) = [1 - (1 - Pdark(ai, a 2 ))(l - r7dct) n ]Pn,+ (ai, a 2 ), (8.7) 

n=0 

with P n>+ (ai, a 2 ) defined by Eq. (3.8) with a + (uo) in place of a(u).] 

As applied to states with photon-number components negligible except when r] dct n «C 1, 
we notice 

1 - (1 - Pdark)(l - 7?dct) n ~ Pdark + (1 - Pdark)^de t n, (8.8) 
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from which we obtain for this case an approximation that simplifies Eq. (8.6) 

oo 

Mi (a) « ^ bdark(a) + (1 - Pd a vk(a))r] dct n]P nt+ (a) 

n=0 

poo 

= Pdark(a) + (1 -Pdark(a))77dct / du a) ' (u) tt(u) . (8.9) 

JO 

C. Detection probabilities 

The most general state involving k modes can be written as 

|^)=J]^ n :POl n (a t 1 ,...,4)|0), (8.10) 

n=0 

where h n is a function of n frequency variables and Pol n (a|, . . . , a\) is a homogeneous poly- 
nomial of degree n in the k creation operators without any annihilation operators. 

Let x stand for any of the modes a,-. The operator for 'detect' for mode x is Mi(x) and the 
operator for 'no-detect' is obtained from Eqs. (8.6) and (3.4) as 

oo 

M (x) = 1 - M 1 (x) = (1 - Pdark (x)) (1 - ^dct) n Pn, + (a:). (8.11) 

n=0 

Any theoretical outcome produced by APD detectors as modeled here is specified by two lists 
of modes, a list J for which APD detectors register 'no-detect,' and a list J 1 for which APD 
detectors register 'detect.' The operator for 'detect' for modes in the list Ji and 'no-detect' for 
modes in the list J is Mi(Ji)M (Jo), where we define 

M (J ) = J] M (x), 

Mxpx) = J] M l {x). (8.12) 
The corresponding probability is then 

Pr(J , Ji) = (VlMiGLjMotJoM- (8-13) 

For expressed in the form of Eq. (8.10) and the APD model that invokes M 1 as defined 
in Eq. (8.6), calculating this probability is surprisingly simple. From Eqs. (8.11) and (3.6) 
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we obtain the effect of M (aj) on a state vector defined by an integral over frequencies of a 
monomial in creation operators 

n n 

M {a) n« f K)|0) = (l-Pdark(a))(l-%ctrn at K')l°) 

3=1 3=1 

n 

= {l-p A ^{a))J[[{l-r ]Aet )a\u j )]\Q,). (8.14) 

3=1 

From this we arrive at the following important 

Proposition: The effect of an operator M (a) for 'no-detection' on a general state is to 
multiply the state by (1 — pdark(a)) and to replace every instance of a creation operator a^(co) 
by (1 — r] dct )a^(ujj). Further this holds if creation operators for modes orthogonal to a enter the 
polynomial Pol„(a|, . . . , a\). 

This generalizes to 'no-detection' of more modes. 

Proposition: The effect of a product of creation operators M (a 1 )M (a2) ■ ■ ■ M (ae) on a 

general state is to multiply the state by 

t 

n^-Pdark^-)), (8.15) 

3=1 

and to replace every instance of a creation operator aj (u), for j = 1, ...,£, by (1 — r] dct (aj)) 
a)(uj). 

For purposes of calculating probability terms (^|M (L)|^), where L is an arbitrary list of 
mode names (as in Sec. 9), we can put this in a symmetric form. Because only like powers 
of creation and annihilation operators appear in terms that contribute to probabilities, we can 
distribute the factor (1 — r] dct (aj)) evenly between the creation and the annihilation operators, 
as follows. 

Proposition: For any list L of mode names, one has 

(V>|M (L)|V) = Y[(l-P<i^(x))(i>'\i>'), (8-16) 

where is the expression obtained from defined in Eq. (8.10) by replacing every instance 
of a creation operator a'j(u), for j = 1, . . . , £, by [1 — ^^(a^^a^Uj), and is obtained 
from (ip\ by the corresponding replacement of a (a;) by [1 — r/d^a^Y^aj^j). 
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Caution: This and the preceding two propositions require that the polynomial in Eq. (8.10) 
contain no annihilation operators. 

For evaluating Mi(Ji), the story is more complicated. What makes Proposition (8.16) work 
is that the substitution of [1 — r7 dct (a : ,)] 1 / 2 a](cjj) for a^{uj) commutes with products, and this 
does not hold for the analogous rule for Mi. From Eqs. (8.14) and (8.11) we find 

n n 

M 1 (a) n^KOlO) = [1 - (1 -p dMk )(l - ?7det) n ] n ot K)l°>- ( 8 - 17 ) 

3=1 3=1 

This proves useful in Sec. 9; however its use is constrained because it cannot be interchanged 
with the taking of products of operators. For this reason we benefit from the following method 
of evaluating Mi(L) in terms of terms of the form M (Lj). Equations (8.11) and (8.12) imply, 
for any set L of mode names, 

Mi(L) = JJ[l-M (x)]. (8.18) 
With this, we express the product of M 1 factors in Eq. (8.12) by 

#(L) 

Mi(L) = ]T(-l)#( x )Mo(X), (8.19) 

XCL 

where for any set S, #(S) denotes the number of elements in S, the sum is over all subsets of 
L, including both L itself and the empty set 0, and we adopt the convention that 

M o (0) = 1. (8.20) 

Example: 

MiK b 2 ) = 1 - Mo(ai) - M (b 2 ) + MqK b 2 ). (8.21) 
Thus we arrive at an equation for the detection operators in Eq. (8.13): 

M (J )Mi(Ji) = (-l) #(X) M (Jo||X), (8.22) 

XcJi 

where J ||X denotes the concatenation of the two lists of mode names. This implies 

(^|Mo(Jo)Mi(JiM = ^(-l) # W(V;|Mo(Jo||XM 

XcJi 

= (-l)#( J o) J2 (-l) #{Jol|X) (V|M (Jo||X)|^). (8.23) 

XcJi 
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D. Effect of time bounds on detection 



In many applications a detector is gated on only briefly. Approximating the turn-on and 
turn-off as perfectly abrupt, we model the effect of gating the detector on for a duration T 
centered at a time t g by use of Eq. (5.9) with the infinite limits replaced by finite times t g — T/2 
and t g + T/2. For instance in Eq. (8.9) we replace J °° dw aj(u)aj(u) by 

rt g +T/2 

/ dta] ± (x,t)a j± (x,t) (8.24) 

Jt g -T/2 

for whichever sense of the ± sign corresponds to propagation toward the detector. Substitution 
from Eq. (5.1) and carrying out the time integration yield 

rt g +T/2 

dt Oj±(x : t)dj±(x, t) 

tg-T/2 



f 

Jo 



°°dco /""dl/ aUu)a j± (^ u ^^-^^ - Sln[(a; ~ U ' )T/2] ■ (8.25) 

JO 7T LO - UJ 1 



[* Look at Mandel [10], first at pp. 573ff; then pp. 69 Iff. Check in Mandel, p. 696 (14.2- 
17), for arguments that dark counts are in principle unavoidable, due to vacuum fluctuations.] 



E. Detection, energy, and photon subspaces 

Let H be the vector space of light states generated by integrals over monomials in creation 
operators acting on the vacuum state |0). Then for a mode a, P n (a)H is the subspace of states 
generated by a weighted integral over n factors, rifc=o a ^ { u k), along with any number of factors 
of creation operators for modes orthogonal to a, all acting on the vacuum state |0). Any vector 
\v) mH can be written as a sum over terms in mutually orthogonal subspaces: 

?1=0 

with |i> n ,a) £ P n (a)H and of unit norm. Then we have (v n , a \v m ,a) = 8m,n- Further, we have 
a(u})\v ,a) = and, for n > 0, a(a;)|i; nj0 ) G P n -i(a)H. It follows that 

(Vw,w') (v\a f (u)a(uj')\v) = ^2(v rija \a t (uj)a(uj')\v nta ). (8.27) 

n=0 
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That is, there are no cross terms. The probability of detection of a mode a, as well as the energy 
in the mode, are often expressed by integrals over products a\u)a{u) or, per Eqs. (8.24) and 
(8.25), aJ(u)a(u'), with the result that both the energy for the mode a and the probability of de- 
tection can be expressed as sums with no cross terms between terms with differing numbers of 
a photons. That is, the expectation energy for energy in the a-mode is h ^ n =o n \C n \ 2 u n , where 
uo n is an averaged frequency that is straightforward to work out. Similarly, using any of the de- 
tection models discussed above, the probability of detecting the state \v) by use of a detector 
involving only the mode a is Pr(Detection of \v)) = J2 n =o | C«| 2 Pr (Detection of \v nA ))- 

F. Preceding the APD detector by a beam-splitter 

The simplest model of beam splitting (which neglects frequency dependence, reflection, and 
mixing of polarizations) expresses an in-mode a in terms of mutually orthogonal out-modes b 
and c by an SU(2) transformation: 

a ] {oo) = (^rans) 172 ^) + (1 - Vtrans ) ( W ) . (8.28) 

Consider an n-photon a-mode in-state defined in Eq. (2.19), 

|^ n ) = (n!)- 1 / 2 / in :at™|0), (8.29) 

where, without loss of generality, h n is symmetric under permutations of its n arguments, 
normalized as in Eq. (2.22). The out-state from the SU(2) transformation expresses the state 
downstream of a beam-splitter as 

Hn) = (n!)" 1 / 2 I U I £ 2 „ s (l - ^a„s) (n - fc)/2 |<M, (8-30) 
k=o \ k J 

where we define the unnormalized state vector 

Hnk) = (h n :b^ n ^)\0). (8.31) 

Now calculate the probability of 'no-detect' for this state by a detector of the b mode, de- 
scribed by the operator M (5) per the preceding APD model: 

Pr (no detect 6) = (V>„|M (&)|^) 
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2 

n 



= ^ransl 1 - ^trans) 



(n-fc) 



k=0 



k 

x (0\(h* n .b k c n - k )M (b)(h n :b^ k S n - k) )\0}. (8.32) 



By Proposition (8.16), this becomes 



1 n in \ 

(Vn|M (6)|^) = (1-Pdark)-^ ^ (1 - ?7trans) <»"*> ( 1 - ^ dct )* 



n! 

fc=0 



k 

x (0|(/i;:6 fe c n - fc )(/i n :6 tfc c t(n - fc) )|0). (8.33) 



From Lemmas (B18), (B17) of Appendix B, along with the normalization of h n we have 

(0\(h* n :b k c n ~ k ){h n :tf k S n - k) )\0) = k\(n - k)\ , (8.34) 

whence we obtain 

(Vn|M (W> = (l-Pdark)^ ( " ] vLns ( 1 " ^ans) ^ ( 1 " V*et )" k\ (n - fc) ! 



fc=0 



fc 



= (l-Pdark)J] " [mrans(l -%ct)] fe (l -mrans)^ ^ 
fc=0 \ / 

= (1 - Pdark)[?7trans(l ~ Vdet) + 1 ~ ^trans]™ 

= (1 - Pdark)(l - ?7trans?7det) n - (8-35) 

This is just the probability of 'no-detect' for the in-state by an a-mode detector, modified by 
replacing rj det by the product r] tTims 7] Act . Thus we have arrived at the 

Proposition: For the APD model described above, we consider a beam splitter for which 
the free input is a vacuum state, and which passes a frequency-independent fraction of en- 
ergy ^trans ; then the probability 'no-detect' downstream of the splitter is obtained from the 
expression for the probability upstream by replacing the detector efficiency r] dct by a reduced 
efficiency i] trans r] det . 

Caution: This proposition holds specifically for the APD model; there are certainly other 
models, such as photon counting, to which it does not apply. 
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9. POLARIZATION-ENTANGLED LIGHT FOR QKD 



Some early models of quantum key distribution assumed bi-photon states [14]. A year ago 
John Schlafer asked how the energy of the light pulse — its multi-photon content — affects the 
statistics of detection. As preparation for Sec. 10, where we start to answer this question, here 
we show some of the possible multi-photon states that have to be considered. 



A. Bi-photon light states 



The general bi-photon light state is given by Eq. (7.5). Here we confine ourselves to a 
special case in which detection statistics are invariant whenever the same polarization transform 
is performed on both the a and b fibers leaving the source. (It will be interesting to check 
experimentally for this invariance.) The polarization transformations are the group SU(2), and 
we refer to bi-SU(2) invariance as invariance of detection probabilities under the following 
frequency-independent transformation, for arbitrary complex u, v such that \u\ 2 + \v\ 2 = 1, 







uai{uj) + va 2 (cu) 








a 2 {u) 




—v*ai(u) + u*a 2 (uo) 




62(a)) 





ubi{uj) + vb 2 {uj) 
—v*bi(u) + u*b 2 (£j) 



(9.1) 



By analyzing the case u — 0, v — 1, one sees that bi-SU(2) invariance requires that h in Eq. 
(7.5) satisfy 

hn(u,u) h 12 (uj,u) 
h 21 (uj,uj) h 22 (u,u) 








l 










-l 






(9.2) 



where (g,g) = 1 but g is an otherwise arbitrary function. Any such state is a single a-photon 
together with a single 6-photon. Sufficiency of this condition to assure bi-photon invariance is 
shown by the substitution defined in Eq. (9.1). 

To explore the dependence of detection probabilities on energy, we need to analyze multi- 
photon states. We confine ourselves to an opening step in this direction by choosing multi- 
photon states that have bi-SU(2) invariance of the form 



\ij) = J2Cn\^n), 



(9.3) 



n=0 
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where 

\^ n ) = f:(a\bt-atb\) n \0) 

/n 
dw x dtii ■ ■ ■ dw n du n f(u u u u . . . , u n , u n ) JJ (a\(uj j)b ] 2 (uj j) - a\(u j)b{(u ^ |0). 

3=1 

(9.4) 

Because the factors under the product all commute with one another, what matters about / is 
the part that is invariant under permutations that swap (luj, luj) with (u k , uo k ). We think of a 
vector ujj = (uj k , Co k ) and denote the average of / over permutations of these vectors by 

S(uj u ...,u n ) f{u u ...,u n ) = ^-y^ /(^tti^tti, • • ■ ,uj nn ,uj nn ). (9.5) 

nl 

This is the first of several symmetries that will be seen to couple polarization entanglement to 
frequency entanglement in a way that affects how detection statistics depend on light energy. 
The terms \ip n ) have the further decomposition 

n 
m=0 

where the unnormalized states \tp nm ) are defined by 

iVVim) = (-l) n ~ m ^ U ^ J dux d£bi ■ ■ ■ du n dw n 

[S(u u . . .,u m )S(u m+1 , . . .,u; n )S(cDi, ...,u n ) /(wi,£>i, . . .,u n ,u n )] 



X 



Vj=l / \j'=m+l 



with the convention on products defined in Eq. (B13), and with the additional symmetry oper- 
ations that act on ujj without acting on ujj, so that S(uj 1: . . . , uo rn ) is defined for m < n by 

S{Ui, . . . , LU m ) f{LUi,CUi, . . . , LU n , LJ n ) 

= f(w n l, d>l, • • • , UJ nrn , UJ m , LU m+ i, & m+ l, ■ ■ ■ , LV n , U) n ). (9.8) 

ml 

Similarly, S{u m+ i, . . . , uo n ) operates on the arguments uij for j — m + 1, . . . , n. 
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Remark: Although S as defined in Eq. (2.20) is not a quantum operator, it is a linear operator 
on a function space, indeed a projection, so that, as an operator, S 2 = S. If h is a function with 
arguments operated on by S then Sh is invariant under the group of permutations over which S 
averages. In particular, Sh is invariant under the swapping of any two of the arguments listed 
in S. That is the first point. The second point is that a transposition of Uj and Uk followed by 
a transposition of ujj and uo k is a transposition of Cuj and Uk- 



B. Effect of a beam splitter 

Generalizing on Eq. (8.28), if we neglect loss, reflections, and unwanted polarization cou- 
plings, we model the effect of a non-polarizing beam splitter (or fiber coupler) that splits say 
mode a>j into modes aji and a j2 by 

a] =ua) l +va) 2 (9.9) 

for some complex u, v such that \u\ 2 + \v | 2 = I. To see the effect of this splitting of a x on the 
state \ip nm ), one makes the substitution for a x in Eq. (9.7) defined by Eq. (9.9). The important 
observation is that no further permutation symmetries enter, so the effect is only to replace 

m 

( 9 - 10 ) 

by 

£ ( m ) (fla{M)) ( ft ■ (9-11) 

k=o \ k J \j=i J \j=k+i J 

C. Effect of polarization rotation 

Unlike beam splitting, the effect of a polarization rotation involves an added permutation 
symmetry, and accounting for this symmetry is essential in calculating the inner products 
needed to arrive at probabilities of detection. This additional symmetry occurs because po- 
larization rotation mixes, for example, an ai-mode with an a 2 -mode. 
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[Part II of this paper contains Section 10, Appendices A through F, and the References, as 
outlined in the Table of Contents.] 
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